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Three-dimensional toroidal precession solitons with a nonzero Hopf index, which uniformly move 
along the anisotropy axis in a uniaxial ferromagnet, have been found. The structure and existence 
region of the solitons have been numerically determined by solving the LandauLifshitz equation. 
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Static and dynamical topological structures with a nonzero Hopf invariant (hopfions) in various models and media 
were discussed in jl], ^, ||, In particular, the magnetization distribution in three-dimensional ferromagnetic materials 
is characterized by an integer Hopf index H. The investigation of the dynamics of the hopfions in these media is not 
only of theoretical interest, but also important for various physical applications, particularly, in view of the prospect 
■ of the creation of fundamentally new memory elements. 

Stable precession solitons with H — (magnon drops) in a uniaxial ferromagnet were found in [j5], However, 
precession uniformly moving toroidal solitons at H ^ were observed only in an isotropic ferromagnet |Q. Their 
stability with respect to perturbations violating the axial symmetry was discussed in The aim of this work is 
to analyze the existence region and structure of toroidal hopfions uniformly moving along the anisotropy axis in a 
c/3 uniaxial ferromagnet. 

The dynamics of the magnetization vector is described by the LandauLifshitz equation; in the case of negligible 
relaxation, it has the form 
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— = - 7 [MxH eff ], H eff = - — , (1) 
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where 7 - is the gyromagnetic ratio (7 > 0). The energy of the ferromagnet is the sum of the exchange energy 

E exch = (a/2) I (d t Mf dr (2) 
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and the energy of the uniaxial magnetic anisotropy 
^ i E ams = 08/2) / (M x 2 + M y 2 ) dr. (3) 
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Many magnetic media are characterized by a high quality factor Q — 0/(4tt) and the contribution of the energy of 
the magnetic dipole interactions is insignificant for them. 

At each point of the K 3 space with a Cartesian coordinate system (x, y, z), the orientation of vector M is specified 
1 by a point on the two-dimensional sphere S 2 in terms of the angular variables and $: 



M x + iM y = MasinQe 1 *, M z = M cosQ. (4) 

Let us consider localized solutions (0), for which 

<S> = ujt + Qp + <l>(r,z-Vt), Q = 6(r,z-Vt), (5) 

where Q S Z, <p - is the polar angle of the cylindrical coordinate system (r, (p, z), and — > at |r| — > 00. Such solutions 
describe three-dimensional precession solitons with a stationary profile, which propagate along the anisotropy axis. 
For simplicity, we consider the configuration of a unit vector field 

n = n(r, z) — (sinOcoscj), sinOsinip, cos9), (6) 

in the coordinate system moving along the z axis (with the velocity V). 

The desired vector field M specifies the mapping M 3 U{oo} — ► § 2 and is characterized by an integer Hopf topological 
index H. If Q ^ 0, the solution given by (||) corresponds to a toroidal hopfion [|[ [To) with the index 

1 />oo />oo 

H = QT, T = — / n - [6Vn x d z n] drdz. (7) 
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Let us represent the energy as a functional of the vector held n: 
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E = I / w E drdz, (8) 
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we = atMo irr 



y2 



(d r n) 2 + (d z n) 2 + ( ^ + ) {n 2 + n 2 



(9) 



In addition to energy, equation (Q) has two integrals of motion: the number of spin deviations (magnons), 
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N = —T I I n z )2nrdrdz (10) 
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and the projection of the momentum of the magnetization held |Tl[ ] on the anisotropy axis, 

Jij poo poo 

P = — / In - [d r n x d z n] Trr 2 drdz. (11) 
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To determine the structure of three-dimensional solitons, we used the same method for minimizing energy func- 
tional (||) with constraint (|l^) as in [T^] , but constraint ( |IT| ) was taken into account through an additional additive 
square penalty function, and an initial held configuration was specified by smooth functions 4> ini and 9 ini , correspond- 
ing to the Hopf bundle with the index H = Q: 

^m = -^ + v + k-, (12) 
2 a 

cos(9 ini ) = 1 - 9 4 ° mh{ ^ e 6(i-~*W), (13) 
2 + cosh(2u) 

where u and v are the toroidal coordinates specified by the relations: 

tanh(u) = (2ar)/(a 2 + r 2 + z 2 ), ^ u < oo, (14) 
cot(v) = (a 2 - r 2 - z 2 )/(2az), < u< 2?r. (15) 

The parameters a and k are determined by the integrals of motion N and P: 



°=^w N - k = h + kF w^' (16) 

where ko ss 3.353 and kp 0.076. 

For comparison, the calculations were performed by the same method for the case of nontopological solitons with 
H = 0, which are stationary ^ and moving ||l3f magnon drops. The initial field configuration in these cases was 
specihcd by the different functions 
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- — + k atan 
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np/Ri, 0^p<i?i, 

7r(2iii - p)/R u Ri< t p<2R l , (18) 
0, p>2R u 



where p = \/r 2 + z 2 , and the parameters k and R\ are expressed in terms of the integrals of motion N and P: 

*-K£*)'(=^)* 

The calculations of the nontopological solitons are in agreement with the data reported in , Jul . 
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FIG. 1: Contours of the angles parameterizing the unit vector n, for the (a) stationary hopfion and (b) hopfion 
moving with the velocity V = 0.134Vb with the same index H — 3 and to — 0.567ujq. The solid and dashed curves 

correspond to 9 = const and cot(<fi) = const, respectively. 



The desired field n(r, z) ensures an extremum of the functional 

J(n(r, z))=E- hwN + VP. 
Using the necessary condition of the extremum, we arrive at a pair of equations 



— J(n(r + \r,z)) 



A=0 



0. 



J(n(r,z + Xz)) 



= 0. 



A=0 



Solving this system with respect to lo and V, we obtain the convenient formulas: 
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V = 



2E, 
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where 



L = aM, 
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(d z n) 2 2irrdrdz. 



(21) 

(22) 
(23) 

(24) 

(25) 



After several thousands of iterations, minimized energy functional (||) reaches a minimum. To test the resulting field 
configurations, we calculated u> and V by equations ( |24| ) and, then, the discrepancy for Landau Lifshitz differential 
equations ([T]). Let us discuss the results. 

Fig.|l| shows the contours of the angles parameterizing vector n. The coordinates are normalized to the characteristic 
length 



(26) 



It is seen that the radius of the central vortex ring corresponding to the value 9 = n, i.e., to the southern pole of 
the § 2 sphere, is larger for the moving soliton. The contours cot((f>) — const are not constructed near small 9 values, 




FIG. 3: Coordinate dependence of the normalized energy density for hopfions of the (a) lower and (b) upper energy 

branches for H = 1, w = 0.567w , V = 0.314V - 
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FIG. 4: Velocity dependence of energy for the hopfion with H = 1 for two values of the precession frequency. 



because this numerical method determining the vectors n does not provide an accurate calculation of the azimuth 
angle when | n z |— *■ 1. However, the shapes of these lines are significantly different for stationary and moving hopfions. 
The localization region of the moving soliton is somewhat larger, but is about Iq as in the stationary case. 

Fig.|| shows the reduced energy e — E/(aM 2 l ) as a function of the ratio V/Vq, where Vq — ^M^y/aft. The plot 
manifests an important revealed property, the existence of two types of moving hopfions with the same H, V and 
oj values at least in a certain velocity interval. Figj^ presents typical distributions of the normalized energy density 
(Te = we / E . A high energy density along the wall of the toroidal surface is characteristic of hopfions of the lower 
energy branch with low velocities, whereas the energy density increases from the wall to the center of the toroid for the 
hopfions of the upper energy branch with high energies. It is also seen in Figs. || and |] that the velocity of topological 
solitons is limited. The limiting velocity of a hopfion decreases with decreasing precession frequency, whereas its 
energy increases. The frequency uj is normalized to the frequency of the homogeneous ferromagnetic resonance: 

uj = 7 M /3. (27) 

The general notions on the structure of the stationary hopfion with H = 1 cannot allow for a comparison of its 
energy with the energy of the corresponding nontopological soliton [jlj], [15| . The features of the structure of the class 
of objects under investigation prevent a numerical analysis of the case with H = 1 and V = jL2j. However, the 
extrapolation of the e{V) dependence to the region of small V values in Fig.|| for H = 1 certainly indicates that the 
energy of the stationary precession topological soliton is higher than that for the nontopological soliton at the same 
frequency u> at least in a certain precession frequency range. 
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